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An analysis of the solutions for the field equations of generalized scalar-tensor theories of gravitation is 
performed through the study of the geometry of the phase space and the stability of the solutions, with special 
interest in the Brans-Dicke model. Particularly, we believe to be possible to find suitable forms of the Brans- 
Dicke parameter uj and potential V of the scalar field, using the dynamical systems approach, in such a way that 
they can be fitted in the present observed scenario of the Universe. 
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I. SCALAR-TENSOR THEORIES OF GRAVITATION 

In a homogeneous and isotropic space, described by the 
Friedmann-Lemaitre-Robertson- Walker metric 



dr' 
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where a is the scale factor and K is the spatial curvature index, 
gravitation can be described by an action of the kind 



S- 



(2) 

where 5™ is the action of usual matter, g is the determinant 
of the metric tensor, a; is a coupling function (which we will 
eventually assume to be a constant, known as Brans-Dicke 
parameter) and V{(j)) is the potential of the scalar field [2]. 
From (2), we obtain for the field equations; 
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where H = a/ais the Hubble parameter and p'" and P'" are 
the energy density and the pressure of the material fluid. 

As usual, we parameterize the equation of state for the fluid 
as P™ = (7 — with 7 a constant chosen to indicate 

a variety of fluids that are predominantly responsible for the 
energy density of the Universe. We can see that through the 
energy conservation equation p'" + 3H{p"^ + P'") = we 
obtain p™ = po/a'^'', with po a constant. 




FIG. 1 : Upper sheet of the of the phase space for a model with lj = 
10 (Brans-Dicke), corresponding to the positive sign in eq. (7) [2]. 
The "hole" in the surface indicates the region forbidden for the orbits 
of the solutions. 



II. THE CASE FOR V = 



AND A' = 



In the referred paper [2], the author proceeds to show the 
phase space allowed for the orbits of solutions for these field 
equations in several cases with different potentials and pa- 
rameters UJ. For example, in the case of vacuum, flat space 
{K = 0) and potential V = ^rn^cp, equation (3) was rear- 
ranged as (making m = 1) 



0, 



which has the solutions 
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The assumption of flat space is required in order to reduce the 
dimensionality of the phase space. 

We want to analyze qualitatively the geometry of the phase 
space {H,4>,(j)), expecting to infer the form of the functions 
u!{(p) and V{(p) to fit better the available data on the structure 
of the Universe. 

The phase space for this situation is composed of a 2-d sur- 
face with two sheets, related to the lower and upper signs 
in eq. (7). Figures 1-3 show the phase space for the choice 
UJ = 10. 

The fixed points for this dynamical system, obtained mak- 
ing H = (f> = Q, are de Sitter solutions, given by Hq = 
±-\/0o/12, with constant Hq and (j)Q. 
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FIG. 2: Lower sheet of the phase space, now corresponding to the 
negative sign in eq. (7). 




FIG. 3; The complete phase space composed of the upper and lower 
sheets linked to each other at the boundary of the forbidden region. 



m. THE CASE FOR V = k(f) AND K = 



write the equation for H: 
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Now, equations (8) and (9) form the system for which the fixed 
points are the solutions Hq = ± ■\/87rpo/30o + A/6 . 

Of special interest is the search for the most adequate func- 
tions tj(0) and V{(j)), that may be more complicated than what 
was assumed until here. 




FIG. 4: Complete phase space for a Brans-Dicke model with a cos- 
mological constant A = 1, energy density po = 2 and constant 
parameter co — 50000, showing two sheets linked by the boundary 
of the forbidden region, as in the precedent case. 



Following other works ([3]-[5]) which give a complete 
analysis of the phase space for Brans-Dicke model with a cos- 
mological constant A (simply making V{4i) ~ Acp in the ac- 
tion), we can illustrate the situation in which uj has a very large 
value and 7 = 0. Therefore, the energy density of the fluid is a 
constant pQ. It should be emphasized that recent observational 
and simulation results seem to favor a scenario very similar to 
this one ([6],[11]-[14]). The solutions in this case are written 
as 



j[<;/)2(A-67i-2) + 167rpo] 



(8) 

With this solutions, we are able to show the phase space for 
a particular choice of constants A, oj and po- 

We proceed to find the dynamical equations system for this 
simple model, as done before. 

Naming A the expression under the root in eq.(8), we can 



IV. CONCLUSIONS 

The method of analyzing the geometry of the phase space 
have proved to be a useful tool in the search for the solutions 
of the field equations of generalized gravity models. Our aim 
is to achieve a complete analysis of the simple model pre- 
sented before (including the stability of the solutions, via Lya- 
punov's direct method [1], in order to investigate further its 
attraction basin) and to apply more sophisticated functions to 
it. 
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